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Abstract 

The phase diagram of a vertex model introduced by P. Di Francesco (Nucl. Phys. B 525, 507 1998) 
representing the configurations of a square lattice which can fold with different bending energies along 
the main axes and the diagonals has been studied by Cluster Variation Method. A very rich structure 
with partially and completely folded phases, different disordered phases and a fiat phase is found. The 
crumpling transition between a disordered and the fiat phase is first-order. The CVM results are confimed 
by the analysis of the ground states and of the two limits where the model reduces to an Ising model. 



*PACS numbers: 05.50.+q (Ising problems); 64.60.-i (General studies of phase transitions); 82. 65. Dp (Thermodynamics of 
surfaces and interfaces). 
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1. Introduction 

Polymerized membranes are two-dimensional networks of molecules with fixed connectivity fl . Their config- 
urational properties and the existence of a crumpling transition between a folded and a flat phase are relevant 
for the behavior of some biological systems j| . Models of polymerized membranes can be defined on regular 
lattices. The constraint of fixed connectivity gives rise to the definition of complicate vertex models which 
offer the advantage that explicit analytical calculations can be performed jj). In these models the bonds 
between the vertices of the network have a fixed length and also, as a further simplification, self-avoidance 
is not considered. In this way the only degrees of freedom of the network are related to its possible states 
of folding and each state can be weighted by a Boltzmann factor depending on the relative angle between 
adjacent plaquettes in the network. 

The model studied in this paper has been introduced by Di Francesco [Q, ^| and concerns the folding 
properties of a square lattice which can be folded along the main axes and along the diagonals. This model 
follows previous studies of the folding properties of a triangular network embedded in d-dimensional lattices. 
In the case d = 2 |^, [?J each triangle of the network can be only in a "up" or "down" state. The entropy of 
this problem is the same as in the three-colouring problem of the honeycomb lattice || . The introduction 
of a bending rigidity induces a first-order crumpling transition which has been studied i n P , |l(| . The model 
with the triangular network embedded in the face-centered-cubic lattice was defined in |lf| and a first-order 
transition has been found also in this case |l2|, 

In this paper we consider the two-dimensional folding of the square-diagonal lattice and study the phase 
diagram of this system in terms of the bending rigidities ks an d Kl relative to short and long edges. Actually, 
the phase behavior of this model was already considered in by means of a transfer matrix analysis, but the 
ground state was not properly investigated, and as a consequence some phases were missing. Here we consider 
the complete set of the ground states of the model and study the phase diagram at finite temperatures by the 
cluster variation method (CVM) |l4[ JlM Eft Jl7[ . This method has already proven to be useful for studying 
vertex models and folding problems |10|, |l2|, |l 3f| . 

The outline of the paper is the following. In the next Section we define the model; the ground states will 
be shown in Section 3. In Section 4 the CVM approximation used in this paper will be described and our 
results will be presented in Section 5. 



2. The square diagonal folding model 

We consider the folding configurations of the square-diagonal lattice where any couple of adjacent triangles 
can be on the same plane or with the two triangles one on the other. Triangles can be folded with respect to 
the main axes of the square lattice or with respect to the diagonals so that in the definition of the bending 
energy different couples of triangles with a short edge or a long edge in common have to be considered. We 
will start the definition of the model by enumerating the folding states in terms of loop configurations. 

2. 1 . Folding and loop configurations 

First fix a reference orientation for all the edges of the lattice in such a way that for each triangle the vectorial 
sum of the oriented edges is zero. There are two choices for this global orientation and one is shown in Fig. 
|I} In any folding configuration each short edge is mapped on one of the four vectors ±ei, ±e2; this mapping 
defines the state of the network. Then consider two triangles with a long edge in common as the pair shaded 
in Fig. 0. This pair can be in two folding states that can be represented in the following way: consider 
another square lattice as in Fig. [| The lines joining the centers of the edges are dashed or full depending on 
whether they are dual to the mapped short edge vectors ±ei or ±e 2 respectively. The obvious observation 
that in each triangle there is a short edge vector in the direction of ±e\ and an edge in the direction of ±e*2 
implies that for each center of the edges of the new square lattice there will be two dashed and two full lines. 
As a consequence the dashed lines and the full lines will form two sets of orthogonal closed loops Q . 

Taking into account also the orientation of e\ and e*2, it comes out that each closed loop has its own 
orientation independently from the others. Therefore the original folding problem is equivalent to a dense 
two-loop problem where a sign representing the orientation has to be attributed to each closed loop. Further 
details can be found in Q. 

Consider now the square plaquettes of this new lattice. In each plaquette there are two dashed and two 
full lines. On each line there can be two possible signs so that there are 2 x 2 4 = 32 states for each plaquette. 
Of course, the sign conservation law along the closed loops implies that only some of the 32 2 configurations of 
two neighboring plaquettes are allowed. The logarithm of the number of all possible states of these plaquette 
configurations gives the entropy of the folding problem. 



2.2. Bending energy 

Here we introduce a bending energy weighting the different folding configurations. We attribute a Boltzmann 
weight e KL or e~ KL respectively to each unfolded or folded long edge. Similar weights e ±KS are associated 
to the two short edge folding states. 

In the plaquette formulation long and short edge bending energies appear in an asymmetric way. Long 
edge bending energy results in an interaction between each couple of adjacent plaquettes in the loop formu- 
lation. An energy — Kl has to be assigned to the network any time a dashed line (or equivalcntly a full line) 
crosses the edge between two plaquettes without bending. In the other case an energy Kl has to be assigned 
to any right angle of dashed lines between neighboring plaquettes (see Fig. ^). 

Differently, short edge bending energy results in different weigths for the 32 basic plaquette configurations. 
In the reference orientation state of Fig. [l], in each square the two short edges on the same diagonal always 
have opposite orientations. These two short edges acquire the same sign for each folding relative to their 
perpendicular diagonal. In this case two triangles are folded and an energy 2k$ has to be assigned to the 
network. In the loop formulation this is equivalent to have a couple of parallel lines in one plaquette with 
the same sign. Similar considerations hold for an unfolded edge and for the other couple of parallel lines in 
a plaquette. Therefore, in total, each plaquette can have an energy equal to — 4ks, 0, or 4k$ @- 

2.3. The vertex representation 

From the discussion above it is clear that our model can be thought as a vertex model defined on a two- 
dimensional square lattice 1? . A vertex variable a x € {1, 32} is associated to each site x of the lattice; to 
each value of a x corresponds one of the 32 vertices introduced above (a possible representation is shown in 
Fig. |). 

Constraints between neighboring vertices are implied by the sign conservation law along the closed loops. 
By using the labelling of Fig. || these constraints can be written in a simple form: we define the two 32 x 32 
matrices jjL h {o-,a') = <5u 2 (<T),u 4 (cr') and p?{p, cr') = Su 3 {u),u-l{o') where a, a' S {1, ...,32} and u(cr) is the vector 
u associated to the vertex a. A configuration a — {<Jx}i? is allowed if /J, h (o~ x ,<T y ) — 1 and /J, v (o~ x ,a y ) = 1 
respectively for any horizontal, (x,y)h, and vertical, (x,y) Vl pair of nearest neighboring sites. Notice that 
on each bond there are only 256 allowed configurations out of the 32 2 = 1024 total ones. We remark that 
when we write (x,y)h,v we suppose the sites x and y lexicographically ordered. 

The hamiltonian of the model is in the form 

H(a) = ]TiTK,)+ E H h (a x ,a y )+ J2 H v (cr x ,a y ) , (2.1) 

zez 2 (x,y)h (x,y)v 

where the two body potential is given by the two 32 x 32 matrices H v (a' , cr) = H h (a 1 a') = 2KLO~v(a),i~v(a') — 
kl for any cr, a' € {1, 32} and the single body potential is given by the vector H s (a) = 4ks if <r = 1, 8, 
H s (a) = —Ks if cr = 25, 32 and H s (a) = otherwise. Finally we can write the partition function of the 
model: 

) J] fi v (a x ,a y ) eM~H(a)} , (2.2) 
where the inverse temperature (3 has been adsorbed in the hamiltonian. 



3. Ground states 

The ground states of the network can be conveniently discussed in terms of loop configurations. First 
consider what happens for ks,kl > at very low temperatures. A large positive kl selects states with 
straight lines so that in a square lattice with periodic boundary conditions there are the dashed and full 
closed loop shown in Fig. ||(a). A positive value of k$ favors alternate sign values on these lines. This 
corresponds to a completely flat state with an energy per plaquette given by E\ = — Aks — 2kl. There are 
8 possible degenerate ground states of this kind and each of them can be realized with two basic plaquette 
states. 

In the sector k$ < 0, kl > elementary triangles want to be folded along diagonals. A positive value of 
kl still select straight lines. All the dashed lines have one sign and the same for the full lines. The network 
is in a state where two triangles sharing a long edge are on the same plane and all the other triangles are 



above one of those two. Also in this case the degeneracy is 8 while one plaquette state is enough for building 
the global state which has an energy E2 = 4k§ — 2kl per plaquette and is reported in Fig. ^(b). 

Loop-configuration ground states for the remaining sectors ns < 0, Kl < and ks > 0, kl < are 
respectively shown in Fig. [|(c) and Fig. |](d). They have an energy equal to E3 = Aks + 2kl and 
E4 = — 4ks + 2kl. They describe a completely folded state with all triangle one above the other, and a state 
with four triangles with a common vertex on the same plane and all other triangles folded above. There are 
needed 2 and 4 different plaquettes for building up ground states in sectors 3 and 4 respectively. 



4. The CVM approximation 

In this section we describe the CVM approximation we used to study the phase diagram of model (2.2). We 
do not enter into all the details related to the CVM technique, we just refer to ju| [IB], lira. 

The CVM is based on the minimization of the free-energy density functional obtained by a truncation of 
the cluster (cumulant) expansion of the corresponding functional appearing in the exact variational formula- 
tion of statistical mechanics. The existence of different horizontal and vertical interactions in the hamiltonian 



(2.1) and the structure of the ground states of the model suggest that one should consider at least a square 
(four-vertex) approximation, that is one should consider a square of four vertices as the largest cluster in the 
expansion of the free-energy functional. Moreover, in order to reproduce the structure of the ground states 
we are forced to partition our lattice into four square lattices with spacing two (see Fig. ||). We denote by 
A, B, C and D these four sublattices and we introduce four square density matrices: Pabcd, Pbadc, Pcdab 
and Pocba- We also need four pair density matrices p A B, Pba, Pen and p DC for the four horizontal differ- 
ent bonds, four pair density matrices p A ci Pca, Pbd and p DB for the four different vertical bonds and four 
single site matrices p A , p B , p c and p D . Finally, we define the sets V = {ABCD, BADC,CDAB, DCBA}, 
C h = {AB. BA, CD, DC}, C° = {AC, CA, BD, DB} and S = {A,B,C,D}, and we write the free-energy 
functional per plaquette as follows: 

f({ Px : X E V}) = \ Exec- Tr* [ Px H h ] + \ Exec* Tr* \p x IF\ + \ Exes ^ \p^'\ 
+1 {J2xev Tr * [Px lo SPx] - Exec- Tr * [Px logp*] 
~ J2xec- Tr A- [Px log/?*] + Exes Tr * \fi* lo SPx}) 
+ J2xev X * (Tr[p*]-1) 

where A x , with X e V , are four Lagrange multipliers ensuring that p x with X S V are correctly normalized 
and for any X G V U C h U C v U S we have denoted by Tr x the sum over all the allowed configurations on 
the set X. 

Following the recipe of the CVM one should find the densities {p x : X £ V{ that minimize the functional 



(4.1) 



(4.1). Hence the next step consists in taking derivatives of the free energy with respect to the square densities. 
In order to do this we must take into account that bond and single site densities can be obtained from square 
densities via a partial tracing. Namely, 

P ab = T^cdPabco Pba — T^dcPbadc Pcd = Tr ' abPcdab Pdc = T^baPocba 

Pac = Tr 'bdPabc d Pbd — Tr ' Ac Pbadc Pca — T^dbPcdab Pdb = ^^caPdcba (4-2) 

Pa = Tr Bj o AB p B = Tr B p BA p c = Tr D p CD p D = Tr D p DC 

Actually, each bond or single site density can be derived via partial tracing of different higher order densities. 
This means that suitable Lagrange multiplier must be introduced to ensure that different partial tracings 
lead to the same result. More precisely, each bond belongs to two different squares (see Fig. |^): we have 
to associate a family of multiplier to each bond to ensure that the same result is obtained by tracing over 
the two plaquettes sharing the bond itself (for instance, we need Tr CD p ABCD = Tr CD p CDAiJ ). We get eight 
different families {X x : X E C U £ v } and each family contains 256 different multipliers, one for each 
allowed bond configuration. Moreover, there exist four different bonds sharing the same single site: to each 
site we associate three different multiplier families {A x>i : X E S, i = 1,2,3}, each of them made of 
32 different multipliers. Hence we have 20 different families of multipliers resulting in a total number of 
8 x 256 + 12 x 32 = 2432 multipliers. 



The functional that we have to minimize is no more the free energy (4.1), but the one in which all the 
Lagrange multipliers are introduced: 



g({ Px : XeV}) = 4f({p x : X G V}) + Tr A [X A:1 Tr B [p BA - p AB \] + 

(4.3) 

"("Tr^B [A AB Tr CD [pcdab ~ Pabcu]] + 

where dots stand for other seven similar bond terms and eleven similar single site terms. Now, let us label 
with a € {1, ...,4608} the allowed 4608 square states. We denote by a x (a) the vertex associated to the site 
X € S corresponding to a. To obtain the equilibrium densities we have to set equal to zero the derivatives of 
the functional (4.3) taken with respect to p Y (a) with Y £ V and a € {1, ...,4608}. In the case Y = ABCD 
we obtain: 

PABco{a) = const x exp {-H h (a A (a), a B (a)) - H v (a A (a),a c (a)) - H s {a A {a)) 
+X AB {a A (a),a B (a)) - A CD (<7 c (a), a D {a)) + X AC (a A {a),a c (a)) 
-X BD (a B {a), <J D (a)) + X AA (a A (a)) + X Ai2 (a A {a)) - X B>1 (a B (a)) (4.4) 
-X Bi3 ((X B (a)) - X ca {a c ,(a)) + X c ^{a c (a))} 

*Pab {<7a (a), u B {a))p AC (a A (a), a D (a))p A (a A (a)) -1 

The equation above is actually a set of 4608 equations. Three similar sets can be found by considering the 
cases X € {BADC, CDAB, DCBA}. The complete set of 4 x 4608 = 18432 equations, together with the 
equations for the multipliers, can be solved by means of the natural iteration method Jl8| , [l9| . The equations 
for the multipliers can be obtained by taking derivatives of (4.3) with respect to the multipliers and are of 
the form 

A x (a x ) = A x (a x ) + constlog^4 ( 4 - 5 ) 

where X E C h U C" U S, a x is an allowed configuration on X and (fi and ip2 are linear combinations of 
higher order densities on clusters Y D X traced over Y \ X. 



5. Phase diagram 

The phase diagram of the model, as predicted by our approximation, has been reported in Fig. ||, where 
open symbols denote second order transitions, and full symbols denote first order transitions. The transition 
lines have been drawn using a limited number of points, since the precise determination of such points is a 
very computationally demanding task. In particular, the transitions between the phase Dl (Disordered 1) 
and the phases F (Folded) and LF (L-Folded) have been simply sketched since it was possible to determine 
their location only with a rough approximation. 

Several phases appear in the diagram. First of all we have four long-range ordered phases corresponding 
to the four possible ground states described in Section ||. Phase Fl (Flat) is stable for large and positive 
ks and kl and represents a flat phase. Phase SF (S-Folded) is stable for large enough absolute values of 
us < and kl > 0. It represents a partially folded phase in which folding occurs mainly along short edges. 
Phase F is stable for large enough absolute values of ks < and kl < 0. It represents a completely folded 
phase. Phase LF is stable for large enough absolute values of ns > and kl < 0. It represents a partially 
folded phase in which folding occurs mainly at long edges. Like the corresponding ground states, all these 
phases have degeneracy 8. 

On the high temperature side of the ordered flat phase Fl we have a small slice of the disordered phase 
D2 (Disordered 2). In the central part of the phase diagram, we have the disordered phase Dl, which has 
larger entropy and larger energy than D2. It is noteworthy that the entropy of this phase at the infinite 
temperature point k$ = kl = is ~ 0.9204 in our approximation, while the estimate in obtained by 
transfer matrix methods, corresponds, with our definitions, to ~ 0.9196. 

Between phases SF and F we have the partially ordered phase PO (Partially Ordered). For large negative 
Kg the transitions between this phase and phases SF and F tend asymptotically to \k^ \ = i In 5+ ^^ ~ 0.412, 



which is exactly the estimate for the critical coupling of the square lattice Ising model in the present 
approximation jljj ]. This is a consequence of a property of the model, which reduces to a square lattice Ising 
model in the limit ks — > — oo. In the loop gas formulation, k$ — > — oo implies that all loops in the same set 
(that is, all dashed loops, or all full loops) must have the same sign. Looking at Fig. || one can verify that 
to satisfy this condition the state of the system must be a mixture of only 2 out of the 32 plaquette states. 
For this pair of allowed states there are four different possibilities (and thus the phase PO has degeneracy 
4): (1,5) (corresponding to + signs on both loop sets), (2,6) (+ on full loops and — on dashed loops), (3,7) 
(— on both loop sets) and (4,8) (— on full loops and + on dashed loops). Hence we have the breaking of the 
loop sign inversion symmetry, which is restored only at the transition to the phase 5. Given a pair of allowed 
states, the hamiltonian reduces to an ordinary Ising hamiltonian, with equal states on adjacent plaquettes 
giving a contribution — kl to the total energy, and different states giving +k^. It is therefore an exact result 
that, in the limit ks — > — oo, there must exist the three phases SF, F and PO (corresponding respectively to 
the ferromagnetic, antiferromagnetic and disordered phases of the Ising model), separated by second order 
phase transitions at |kl c | = | In (l + \/2) ~ 0.441. 

A similar situation occurs in the limit kl — > — oo. Looking at Fig. |]one sees that in this limit both the 
full and the dashed lines form the smallest possible square loops. It follows that the model reduces to two 
decoupled Ising models, where the Ising variables are the loop signs and two loops (both full, or both dashed) 
having parallel edges in the same plaquette interact with an energy ±2k$- Therefore the phases F and LF 
(corresponding respectively to the antiferromagnetic and ferromagnetic phases of the limiting Ising model) 
will undergo second order phase transitions towards the disordered phase Dl at \ksc\ = j m (l + v2) — 
0.220. The estimate for this value in our approximation is of course \kl\ = \ In 5+ /^ — 0.206. 

Comparing our phase diagram with that by Di Francesco S one can see several striking differences. 
First of all, among our ordered phases, only the completely flat and completely folded ones were reported 
in ||, and no low temperature transition was found among them. In addition, no intermediate phase (like 
our phases PO and D2) was found between the ordered and the disordered phases. We note that the CVM 
results are fully confirmed by the ground states analysis and the limits k$ — > — oo and kl — > — oo. 

In conclusions, we have shown that the CVM techniques can be adapted to study very complex vertex 
models as the one considered in this paper representing the folding configurations of a square-diagonal 
lattice. As in the case of the folding of the triangular lattice, a first-order crumpling transition between 
the flat phase and a disordered phase has been found. The extension of this analysis to the d-dimensional 
folding problem appears not easy from a numerical point of view. 
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Fig. 1: The square-diagonal lattice. 
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Fig. 2: The mapping onto the vertex model. 
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Fig. 3: A possible representation of the vertex variables. Each vertex is labelled by a variable 
G {1, 32}. Following [0 a vector u — (ui, U2, 113, 114) G {0, 1, 2, 3} 4 and a scalar v £ {0, 1} are 
associated to any a £ {1, . . . , 32} (see 0] for more details). 
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Fig. 4: The ground states depicted in figures (a), (b), (c) and (dlcorrespond, respectively, to the 

phases Fl, SF, F and LF in Fig. g. 
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Fig. 5: Partition of the lattice used in the CVM calculation. The eight bonds and the four 

plaquettes are explicitely depicted. 




Fig. 6: The phase diagram: open (full) symbols denote second (first) order transition points 
calculated via CVM. Solid lines are eye-guides, dashed lines are suggested by the limits ks, 

Kl — > CO. 



